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Suspensions of rodlike objects in a liquid are encountered in many areas of science and technology 
and the need to orientate them is extremely important to enhance or inhibit certain chemical 
reactions between them, other chemicals or with the walls of vessels holding the flowing suspension. 
Orientation control is feasible by altering velocity and nature of the liquid, its flow and the geometry 
■ of the channel containing the flow. In this work we consider the simplest possibility of orientation 

control with flow in two dimensions on the basis of the orientation distribution of rodlike objects. 
The simple differential equation satisfied by the orientation probability density function is derived 
£SJ ■ and its solution discussed. In addition, we show how birefringence and dichroism ("Maxwell effect") 

^ ' of the suspension provide a direct experimental test of orientation control. 
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In many physical, chemical, biological processes, the behavior and orientation of rodlike objects (such as fibers, 
nanotubes, nanowires, DNA, macromolecules) suspended in a flowing liquid affect the transport, rheology, chemical 
and hydrodynamic characteristics of the suspension. Orientation control for the sake of aligning or separating the 
rodlike objects (RLO), enhancing reaction between them, with other chemicals or vessel walls is important in many 
areas of science and technology such as sedimentation, blood flow, pulp and paper, polymer processing, liquid crystal 
flow, microfluidic devices, ferrofluids... 

This work is relevant to students who have completed an undergraduate Statistical Physics course of the Reif [l[ 
fS |' level and are interested in applications of stochastic processes or graduate students whose level corresponds to Landau 
and Lifshitz course [|[ and are interested in the general applications of Fluid Mechanics. The history of dilute solutions 
| of RLO (from the polymeric point of view) is detailed in Chapter 8 of Doi and Edwards book " The Theory of Polymer 
Dynamics" 4] as well as in a review paper by Forest et al. Q. Students interested in various other materials and 
' applications are encouraged to look at the "Soft Matter" four volume series edited by G. Gompper and M. Schick 
Q\ , In this work, we are interested in dilute suspensions of rigid and long RLO with a negligible cross-sectional area and 
concentrate on the study of orientation control by flow in two dimensions. Our goal is to investigate in the simplest, 
yet physical case, and with a self-contained pedagogical point of view, the effect of liquid shear flow on the rotational 

■ diffusion of a dilute suspension of non-interacting RLO in a planar contraction. We derive in a straightforward manner 
the main equation describing the statistical distribution of orientation in two dimensions (2D), solve it analytically 

0^ . and show how orientation control might be detected optically with a measurement such as birefringence (double 
refraction) and dichroism of the suspension. 
, Objects suspended in flow undergo two types of motion: smooth motion due to the average fluid velocity field 
and erratic random motion produced by the fluctuating fluid velocity, temperature and inertia driven motion. The 
resulting change in the suspension microstructure can have a significant effect on the mechanical, thermal, optical, 

■ electrical, magnetic and chemical properties. The equation that accurately models the orientation state of non- 
interacting RLO in hydrodynamic nonhomogenous flow depends on the rotational Peclet number a that represents 
the interplay between the randomizing effect of temperature induced rotational diffusion and the orientating effect of 
streamwise mean rate of strain due to the liquid flow. 

Two types of flow are of particular importance: simple shear flow and elongational (extensional) flow. Simple 
shear flow corresponds to a velocity profile where the gradient in the fluid flow velocity is constant. In the case of 
elongational flow the sample is compressed in one direction and elongated in the other. In this work we consider 
a dilute concentration of non-interacting RLO suspended in a liquid performing shear flow exclusively in a planar 
contraction of width L (see fig. [J). 

The rotational Peclet number a measuring the relative strength of hydrodynamic interactions and Brownian forces 
is defined mathematically as the ratio of hydrodynamic shear flow and the rotational diffusion constant: 
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FIG. 1: 2D geometry of the shear flow and the RLO of length I making an angle with the flow direction along x. The shear 
rate of the flow is 7 = • The fluid velocity profile v x is shown explicitly as a function of the coordinate y. In the frame 
moving with the rod, the velocity field is antisymmetric (resp. symmetric) above (and resp. below) the rod. 



a = w r w 

7 is the (hydrodynamic) shear rate and D r the (thermal) rotational diffusion coefficient (see the Appendix for an 
explanation about the difference between ordinary diffusion 111] and rotational diffusion). 

The flow shear rate defined by 7 = with v x the velocity field in the x direction of the flow (see fig. [1} is related 
to translational degrees of freedom whereas D r , the Brownian diffusion coefficient, governs the rotational motion of 
the RLO around its center of mass and hence relates to rotational degrees of freedom (see Appendix). 

The RLO suspended in the liquid are under the action of two orientating forces: one that comes from the hy- 
drodynamic flow and the other due to thermal effects. Despite the apparent simplicity of the problem, its physics 
is extremely rich and one might stumble upon some unusually complex behaviour; for instance the following set of 
states may be observed: 

1. An aligned state where the RLO is aligned with the flow at a fixed angle; 

2. A tumbling state, in which the RLO lies in the shear (x — y) plane and rotates about the z axis. 

3. A wagging state, in which the RLO lies in the shear plane, but oscillates between two values. 

4. A Kayak-tumbling state, equivalent to the tumbling state, but in which the RLO is out of the shear plane. 

5. A Kayak- wagging state where the RLO is out of plane, but the projection of the RLO on the shear plane 
oscillates between two values. 

6. A mixed state, typically found close to the boundaries between wagging and tumbling states and where the RLO 
undergoes both oscillation and complete rotation. 

7. A Chaotic state, in which the RLO angle with the flow behaves with time in an apparently random way. 

While the above describe the time dependence of the RLO orientation, it would be rather interesting to evaluate the 
statistical distribution of the angles the RLO makes with the flow. Specializing to 2D, Boeder Q is the first to have 
studied the RLO stationary angular distribution problem in a flowing liquid from a theoretical and experimental point 
of view. In the shear flow case, he derived an ordinary differential equation (ODE) satisfied by P(9), the stationary 
PDF (see Appendix) describing the average orientations of the RLO assuming that the motion of RLO, of negligible 
cross-sectional area occurs in the plane of the flow, without any boundary conditions. 

The resulting ODE reads: 

^P(9) + ^[asin 2 (9)P(9)]=0 (2) 

In order to derive the above ODE, the RLO angular time evolution induced by shear flow is needed. We assume 
that the effects of gravity and inertia on the RLO are negligible as well as any disturbance motion resulting from 
the RLO presence. Taking a unit vector u along the RLO axis, u = xcos9 + ysm9, its time derivative is given by 
<it = ( — ^ sm ^ + ycos9). This yields the angular time evolution of the RLO as z4| = u x 4| where z is the unit 
vector perpendicular to the x — y plane. In addition, the fluid velocity field v = v x x is along the x direction and 
depends solely on the coordinate y. The above simplifying assumptions imply that the RLO follows the surrounding 
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fluid velocity gradient and that the time derivative of u along the x direction (see fig. [J) is proportional to the shear 
rate (v x change along y) and the u y component (see additional note 3). This means ^jf- — ju y , i.e. % = jsinOx. 
Therefore, the RLO angular time evolution induced by shear flow becomes: 



(3) 

Thermal disturbances to the RLO motion, are accounted for by adding to the above a random term proportional 
to a time-dependent White noise amplitude £(t) possessing the following statistical properties: 



< >= and < f (*)£(*') >= 6{t - t') (4) 

where < £ > is the statistical average of £, meaning that £(t) has zero mean and zero correlation with another 
amplitude occuring at a different time. The angular time evolution becomes: 



d9 = Adt + VB£(t)dt (5) 

with A = — 7sin 2 (#) and B a constant independent of time and the angle 9. This implies that the deterministic 
hydrodynamic forces tend to act on the RLO rotating it in the shear flow with an average angular speed 4£ while the 
fluctuating term originating from thermal effects produce random disorientations. 

The differential equation Adt + y/~B£(t)dt in the variable 9 containing the random term £(t)dt is called a Langevin 
equation [l[ that can easily be transformed into a ID Fokker-Planck equation (see the derivation in Chapter 15 
of Reif 's book [l[ or the paper titled " Simplified Derivation of the Fokker-Planck equation" by Siegman Q in this 
Journal) : 



Straightforward comparison with the ODE given by eq. [5] leads to B = 2D r and in the limit of no shear (a = 0) we 
recover the stationary planar rotational diffusion equation D r ^P{9) = (see Appendix). 

Having derived the Boeder ODE from simple assumptions, this work is meant to provide a solution in closed form 
for the ODE (eq. ^ using standard analytical methods (pertaining to Freshman calculus level) to obtain P{9) for any 
Peclet number a and assess later on the impact with regard to orientation control from the measurement of optical 
properties. 

This paper is organised as follows: In Section II, we present an exact analysis of the ODE (see eq. ^ to obtain the 
PDF, for a wide range of a £ [10~ 4 — 10 s ] . In section III optical properties are evaluated to test liquid flow effect 
on orientation control. Conclusions are given in Section IV. The Appendix explains in detail the difference between 
ordinary and rotational diffusion. 



II. ANALYTICAL SOLUTION 



In this section, an analytic procedure for deriving the solution of the ODE (see eq. \2§ as the probability density 
function (PDF), P{9) is presented. The PDF solution of the ODE is 7r-periodic since the RLO's are indistinguishable 
when orientated at 9 or 9 + %. Since the PDF is 7r-periodic one might write, in principle, a Fourier series solution 
valid for small and large values of a. 

Since we are interested in finding an exact solution we differ from this approximate approach by viewing the problem 
as a boundary condition one, more precisely: 



P(0) = P(tt) and P'(0) = P'(tt) (7) 



where P'{9) = 

Since the PDF ought to be normalised over the interval [0, 7r]: 



f P{9)d9 = 1 
Jo 



(8) 
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Consequently, the determination of the PDF is a constrained periodic boundary value problem eq. [2l the constraint 
originating from the normalization condition. 

In order to get closed form solutions, we write eq. [2] as a first-order ODE: 



P'{6) + asm 2 (9)P(9) = C (9) 

with initial condition: P{6 = 0) = P(0). This is mathematically sound provided the initial value P(0) and the 
constant C are known for any value of a. 

The constant C appears in eq. [5] and is equal to P'(0) by substituting 9 — and assuming finiteness of P(9 — > 0) 
in eq.E). 

The formal solution of eq. [5] may be given generally as: 

, „ „ r a,sin(20) ,, f . a,sin(2:r) , 
P(0) = Cexp[-(— ^-0)]^ exp[--(— ^ - x)]dx (10) 

The lower limit — oo is not obvious since the problem is defined over the angular interval [0, w] . Analytically, the lower 
limit — oo is the only possibility compatible with the boundary conditions given by eq. [7] since the factor exp(aa;/2) 
appearing in the integrand as given by eq. [TU] enlarges the angular interval from [0, tt] to the interval ] — oo, 9]. The 
reason behind this transformation comes from imposing a zero integration constant that can be obtained only when 
x — > — oo given the factor exp(aa;/2) and the positive sign of a. 

Performing a change of variables, the solution may then be written as: 

P{0) = ( — ) exp[^ sin(20)l f°° exp(-x) exp[ £ sin(— - 29)}dx (11) 
V a ) 4 J 4 a 

The form in eq. [11] is superior to the previous form given by eq. 1101 since the the exponential factors exp(#) and 
exp(— 9) terms (that can vary rapidly by several orders of magnitude) are not used (Toj |. The constant C is determined 
from the normalisation condition of the PDF (eq. [5]) yielding: 



1 2 f~ f°° a 2x 2x 

— = - I d9 dx exp(-x) exp[(-)sin( — ) cos(26» )] (12) 

C a J_k J 2 a a 

From the general expression eq. 1111 the initial value P(0) is given by: 



/2<7\ I" 00 n At 

P(0) = I — / exp(-x) exp[ - sin( — )]dx (13) 
\ a J J Q 4 a 

The PDF depends on a which we may want to vary over several orders of magnitude. 

Surprisingly, the difficulty in solving the ODE stems from the fact that its nature may be modified when a increases, 
with the constraints of 7r— periodicity and PDF normalization maintained. For small values of a the PDF is flat as 
observed in fig. [2] and this can be shown with a simple argument. 

When a approaches zero, we should satisfy the pure rotational diffusion equation in the plane given by D r 4s^P(0) = 
whose solution is P(9) — a9 + b with a, b constant and determined from 7r— periodicity eq. and normalisation 
condition eqJS] yielding a = 0, b = —, 

In the opposite case of large a (as in fig. [3]) the PDF displays a strong peak around an angle 9 max indicating that 
orientation of the RLO is setting in. 



III. OPTICAL PROPERTIES 



In general, even if a solution of RLO's is optically isotropic at equilibrium, under flow conditions, the distribution 
of orientation will be affected in a way such that the optical properties become anisotropic because of the existence 
of some privileged direction due to the flow. This phenomenon is called flow birefringence or Maxwell effect [4|. 

Shear-induced birefringence and dichroism measurements of suspensions of RLO consists of detecting phase and 
intensity variation (due to the flow) of light travelling across the sample. 
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FIG. 2: Analytical PDF as a function of 6 for small a=0.01, 0.1. and 1.0 normalised over the interval [— tt/2, tt/2]. The isotropic 
(fiat) angular distribution whose value is •=■ over the interval [—tt/2,tt/2] is reached as a decreases toward zero. On the other 
hand, when a increases, a wavy structure develops about the flat distribution with a maximum appearing around 7r/4. 
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FIG. 3: Analytical PDF as a function of 6 for large q=10.0, 100. 
increases, the PDF becomes sharply peaked around an angle 6 ma:l 



and 1000. normalised over the interval [—tt/2, tt/2]. 
that approaches zero. 



As a 



The flow produces a modification in the real and imaginary part of the optical index An', An" of the suspension. 
An' is extracted from the phase difference between the transmitted and reference signal whereas An" is extracted 
from the attenuation of light intensity. An', An" are both proportional Q to a single quantity, S the nematic order 
parameter encountered in the physics of liquid crystals. 

S is for the unorientated (disordered) state and 1 for the fully orientated (ordered) state. Writing S = a < 
cos 2 6 > +b (with a, 6 constants to be determined) and noting that in the ordered state < cos 2 9 >= 1 since 9^0 
whereas in the disordered state < cos 2 9 > = 1/2. Hence we find S = 2 < cos 2 9 > — 1 =< cos 29 > [121 ]. 

When S =< cos 29 > saturates to one, perfect alignment is reached i.e. complete orientation control is achieved as 
observed in Fig. 2] for very large values of the Peclet number a. 

When S saturates to 1, the orientational PDF becomes sharply peaked around 9 max the value of the angle that 
corresponds to the PDF maximum as the Peclet number increases. At the value 9 — 9 max , the first-order ODE (eq. [9]) 
writes: 

a sm 2 (9 max )P(9 max ) =C (14) 
Since, it is required that the PDF should be normalised according to eq.[5]when a is large, we get after using eq.fHl 



9 m axP(9 max ) ~ constant (15) 

The above is the area under the PDF curve since the distribution function becoming sharply peaked around 9 max 
when a is large, leading to a PDF approximately triangular in shape with a height P(9 max ) and a base equal to 
29 max - This implies that 9 max is close to the standard deviation of the PDF P(9) when a is large. 
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Using eq. [14] and eq. [T5J we find the following leading behaviour: 



' 1/3 , P(e max )~a 1 



/3 



(16) 



Hence S behaves for large a as: 



5 = / P{9)cos29d9 



P{9) cos29d9 



(17) 



which yields: 



S rs -' QmaxPifimax) COS 2 $77 



(18) 



implying that S 1 behaves as cos(20 ?i 



2 cos 2 



found in fig. |4] and obtained from full numerical integration. 



5 ) — 1 ~ 1 — 2a 3 which is close to the asymptotic behaviour 




FIG. 4: Nematic order parameter S =< cos(2#) > versus Peclet number a spanning the interval [10 4 — 10 8 ]. For small values 
of a we have a "disordered" state characterized by S ~ whereas 5* approaches 1 ("orientated" state) for large values of a. 

— - 

The asymptotic regime 1 — 2a 3 is shown for comparison. 



Since 9 max is a measure of the standard deviation (for large values of a) in the fluctuations of the angle 9 about 
the main flow direction, the RLO orientation controllability with flow can be analysed with the variation of the angle 
dmax with the Peclet number a. 

In fig. [5] the variation of 9 max with the Peclet number a is depicted. From the figure, one infers that orientation 
controllability is achieved when 9 ma x drops below a small specified value. For instance, one has to have a ~ 400 in 
order to achieve 9 max < 0.1 radian (see fig.0). 



IV. CONCLUSIONS 



This work is the simplest illustration of the interplay between shear (due to liquid flow) and rotational diffusion 
(due to thermal agitation) on the planar orientation of a dilute concentration of rigid non-interacting thin and long 
RLO suspended in a liquid shear flow in a simple geometry. 

The PDF describing the RLO orientations is analytically determined in 2D and accurately evaluated as a function 
of cv, the rotational Peclet number given by the ratio of shear rate and rotational diffusion constant. 

The solutions found for the PDF in the bulk of the flowing liquid for arbitrary values of a confirm the validity of 
the approach used for small values of a. In this regime flat angular distributions (PDF) are expected: this means we 
have isotropic PDF due to rotational diffusion dominated motion. 

On the other, when a increases, we expect sharply peaked angular distributions (PDF) the peak occuring around 
an orientation angle 9 ma x- this means we are in the shear dominated motion case and the PDF is anisotropic with a 
single mode. 
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FIG. 5: Variation of the angle 8 m ax (in radians) with the Peclet number a. For small values of a, 9 max saturates at a value 
of about 0.8 radian that is close to j as seen in fig. [2] When a is increased beyond 400, 9 ma x decreases as a~ 1//3 dropping to 
small values < 0.1 radian indicating that orientation control has started to set in. 



The impact of a on orientation control and its signature through optical properties is examined and the variation 
of the nematic order parameter S is obtained for a wide range of values of a spanning the interval [1CT 4 - 10 8 ]. 

Since Boeder's pioneering work, many improvements have been made to remove restrictions on the cross-sectional 
area and size of the RLO, introduce internal vibrational and rotational degrees of freedom within the RLO Q, consider 
rotational diffusion in 3D, turbulent flow [13], more complex flow geometries [4J and interactions between the RLO's. 

The Boeder equation provides orientation control requiring a wide variation of the Peclet number whereas practically 
a much smaller range is demanded. From fig. [5] the transition from 9 ma x ~ 7r/4 (almost flat and isotropic behaviour 
leading to no alignment) to the case 9 max ~ (peaked behaviour, anisotropic meaning controlled orientation) spans 
an interval a of [0.1 — 400.] in order to switch from no alignment to orientation control in the 0.1 radian range. If, on 
the other hand, we require 9 max < 0.01 radian, then q~6x 10 5 which is a very large value. 

Practically speaking, a small interval of variation of a is required in order to achieve a fast and precise control of 
alignment (with a small standard deviation), hence the question: 

Would it be possible, within a model similar to this one and under certain conditions to vary slightly a (by a few 
units only) and achieve a standard deviation smaller than 0.01 radian? 

Another question is: Would it be possible, within a model similar to this one and under certain conditions to vary a 
and obtain several peaks (multi-modal) behaviour with several orientations, in other words make a transition with a 
from single-mode to multi-mode behaviour? 



V. APPENDIX: ROTATIONAL DIFFUSION 



Ordinary diffusion describing the random walk of a particle in space is given by: 

dffit) 



dt 



= DAf(r,t) (19) 



with f(r,t) the probability function for the particle to be at position r at time t and A/ the full Laplacian of /. 
D is the ordinary (translational) diffusion coefficient. In rotational diffusion, a particle wanders on the surface of the 
unit sphere such that the diffusion equation becomes: 



2IM. = DrAemt ) (20) 

The angular Laplacian Ag and the probability function / = f(0,t) depends now solely on the angular coordinates 
represented collectively by 0. Ag contains (second-order) derivatives with respect to 9 with the radial variable fixed 
(r = 1) and with the rotational diffusion coefficient D r replacing the ordinary diffusion coefficient D. 

In the 2D case, the rotational diffusion equation on the unit circle reads: 



df d 2 f 
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The stationary PDF P{&) studied in this work is obtained from the long-time limit: P{G) = lim^oo f(9,t). 

In contrast to translational diffusion, D r units are not Length 2 /Time but simply 1/Time since the random walks 
on the unit sphere are purely angular with no dependence on any characteristic length. Similarly to translational 
diffusion where the average displacement squared < [r(t) — r(0)] >^ Dt at time t, in rotational diffusion, we have 
< [u(t) — u(0)] 2 >~ D r t with u(t) representing the angular position 6{t) at time t (while radial position r = 1 is kept 
fixed for all times). 
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